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ABSTRACT 

We compute the effective action in terms of the Polyakov loop for the 3- 
dimensional pure fundamental- adjoint SU{2) lattice gauge theory at non- 
zero temperatures using the strong coupling expansion. In the extended 
coupling plane we show the existence of a tricritical point where the nature 
of the deconfinement transition undergoes a change from second to first or- 
der. The resulting phase structure is in excellent agreement with the Monte 
Carlo results both in the fundamental and adjoint directions. The possible 
consequences of our results on universality are discussed. 
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1. INTRODUCTION 



The deconfinement transition in gauge theories is well known to be inti- 
mately connected to the spontaneous symmetry breaking of the correspond- 
ing center symmetry with the Polyakov loop (L) as the order parameter. As 
a consequence of universality, various gauge models in [d + 1) dimensional 
Euclidean space get related to much simpler (d) dimensional scalar or spin 
models near their critical fixed point [|I], ^. The confining and deconfining 
phases in the language of spin models correspond to disordered (< L >= 0) 
and ordered phases (< L >^ 0), respectively. The effective action written 
in terms of the Polyakov loop, after integrating out all the other degrees of 
freedom, contains all the relevant information about the phase structure of 
the theory. These effective actions for SU(N) gauge theories and their uni- 
versality with spin models have been well studied in the past using analytical 
tools like large N 0, |, || , mean field |^, ^ , strong coupling |||, |^ and weak 
coupling approximations as well as by Monte Carlo techniques [^, pT[1 . 

On the other hand, regulating the theory on a lattice, it is expected that 
the various choices of gauge invariant lattice actions, which in the naive clas- 
sical continuum limit (i.e ^ — 0, a — 0) reduce to the Yang Mills theory 
also correspond to working within the same universality class. However, in 
the absence of full non-perturbative renormalization group equations, differ- 
ent lattice models expected to be in the same universality class have been 
explicitly tested for universality and their scaling behaviours explored both 
analytically and by Monte Carlo simulations. A class of lattice models, with 
a rich phase structure and therefore well studied in the literature is the pure 
fundamental-adjoint SU {N) Wilson action defined by the partition function 



Here F, A denote fundamental and adjoint representations, respectively, 
and Up is the ordered product of the four directed link variables U^{n){& 
SU{N)) which form an elementary plaquette. The sum over P denotes the 
sum over all independent plaquettes of the lattice. / Jl/^.n dU^{n) is the prod- 
uct of group invariant Haar measure over all the links of the lattice. This 
model represents a particular choice of lattice action from a more general 




(1) 
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class defined by summing the above single plaquette action over all the rep- 
resentations of the gauge group with corresponding couplings. As in the 
naive classical continuum limit, the theory is only sensitive to the Lie alge- 
bra, thus the above class is believed to be in the same universality class as 
that of Wilson action. In fact, the universality between (|l|) and the Wilson 
action [i.e, (3a = 0) has been explicitly demonstrated in perturbation theory 
near the critical fixed point and, as expected, the two differ only in their 



scale factors (A). 

In the special case of SU(2), the phase structure and the scaling properties 
of the above model have been extensively studied by Monte Carlo simulations 
[ll, |14|, |15| as well as by analytical methods [||, |16|, [l^, m. Along the 



Pf = axis, it describes the SO (3) gauge theory with a first order phase 
transition at ~ 2.5. At = oo, it corresponds to the Z2 lattice gauge 
theory with a first order phase transition at = | ln(l + ^/2) ~ 0.44. From 
the /3f = and Pa = 00 axes, the above two bulk transitions extend into 
the [jSp, Pa) plane, meet at a point and then continue as another line of bulk 
first order phase transitions. The fact that the latter line ends at a point 
in the phase diagram allows one to bypass all the bulk singularities without 
losing confinement in SU(2) gauge theory. This model therefore retains its 
strong couphng confining property also in the continuum and is expected to 
be in the same universality class as that of the Wilson action {Pa = 0). 
However, as the above phase diagram was mainly established using Monte 



Carlo simulations on relatively small lattices, we recently argued ||2y, ^1 



22 1 that it is incomplete and besides the bulk transitions discussed above, 
one must also worry about deconfinement phase boundary in the extended 
coupling {Pf,Pa) plane. Monitoring the Polyakov loop expectation values 
< |L| > and its susceptibility by Monte Carlo techniques, it was observed 
that the known second order deconfinement transition at Pa = moves into 
the phase diagram as Pa is increased and eventually joins the first order 
bulk transition line at its endpoint around Pa = 1-0. The above merging 
together of the deconfinement transition line with the previously claimed bulk 
transition line lead to a paradox as their nature as well as scaling behaviours 
are completely different. Even more surprisingly, it was observed that 

1. For /3a < 1, the model was in the universality class of the 3-dimensional 
Ising model as predicted by the universality hypothesis ^ How- 
ever, the deconfinement transition was found to be discontinuous above 
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/3a ~ 1 showing a definite qualitative change in the properties of the 
theory and the existence of a tricritical point. 



2. Absence of the previously claimed first order bulk transition. There 
was good evidence that the above transition was misidentified in the 
Monte Carlo simulations and is, in fact, first order deconfinement tran- 
sition mentioned above. It was found that the first order portion of 
the deconfinement transition lies exactly along the previously known 
bulk transition line and there was no evidence of two transitions. The 
plaquette susceptibility does show a peak along the above transition 
on smaller lattices, however, on going to larger lattices, this peak was 
found to decrease considerably with the total 4-volume of the lattice. 



To understand the above unexpected and puzzling lattice Monte-Carlo 
results and their possible implications to the continuum physics and uni- 
versality, one needs to study and derive these results analytically. Such a 
study can yield further insight and a global understanding of such phenom- 
ena for SU (N) lattice gauge theories in general. Analytical results also have 
the advantage of being free from finite size artifacts. With this motivation, 
in this paper we study the deconfinement transition and its nature in the 
above extended SU{2) model in the strong coupling expansion. We explic- 
itly compute the Landau Ginzburg effective action in terms of the Polyakov 
loop and show the existence of the tricritical point in the extended {(3p,i3a) 
plane where the the deconfinement transition changes its order (Fig. 2). 
Moreover, the above effective potential reproduces the Monte Carlo phase 



diagram of pO|, ^ with good accuracy (Table 1). It is also found that 
within the strong coupling regime the tricritical point moves in the direc- 
tion of {Pf = oo, /3a = oo) as A^^ increases. This possibly indicates that 
the continuum limit of the theory may not be affected by the existence of 
the tricritical point. Surprisingly, we find that the predictions of the critical 
couplings along the Wilson line (3a = are quite close to their corresponding 
Monte Carlo values (Table 2), even well inside the scaling region established 
by numerical simulations p4|. In the above case of pure Wilson action, the 
deconfinement transition has also been studied by Polonyi and Szlachanyi 
[§] using strong coupling expansion. However, our techiniques differ from 
theirs and our results in this special case {(3a = 0) are different and closer to 
the Monte-Carlo findings. Following the work [EO, pT|, a more recent paper 
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2^, studied the phase diagram of the fundamental-adjoint model for SU{3) 
lattice gauge theory at nonzero temperature by Monte Carlo simulations. 
They reported results consistent with the usual universahty picture. We will 
briefly summarise and discuss their results at the end. 



1. THE LANDAU GINZBURG MODEL 



We want to compute the partition function in terms of the Polyakov loop 
after integrating out all the spatial link degrees of freedom in (H). As usual, 
we impose periodic boundary conditions in the temperature direction and 
assume its length to be Nr. Defining W{n) = Y[no=i Uo{no, n), the Polyakov 
loop can be written as, 



L{n) = ^TrU^ (n) (2) 

The partition function (|l|) can be trivially written as product over the 
plaquettes. 



n JdU.in)] n (l + ) (l + Gp^) > 



(3) 



.Links 



Plaquettes(p) 



where. 



E 



and 



m„=l ""P 



oo 1 

F — 

^ n 1 



Pa n2 



(4) 



(5) 



Here Sp denotes the standard one plaquette action for the gauge group 
SU{2). 

The computation of the effective action in terms of the Polyakov loops 
get considerably simplified by noticing that after integration over the spatial 
links in an arbitrary cluster of plaquttes in (H), 
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1. Due to gauge invariance of the partition function and the periodic 
boundary conditions, the non-trivial dependence on Uo{no, n) will come 
only if the cluster contains at least one set of plaquettes forming closed 
loops in the temporal direction. Moreover, the above dependence can 
come only implicitly through the Polyakov loops[j] 

2. The rest of the cluster's contribution to the partition function is ei- 
ther zero or proportional to some powers of fundamental or adjoint 
couplings. 

Therefore, after the above spatial link integrals are performed the Landau 
Ginzburg effective action Sejf can be defined as, 

Z = J]Jdfi{L{n))exp-Seff{L). (6) 

n 

Here dfi{L{n)) is the effective measure in terms of the Polyakov loop and 
is computed below. For matter of convenience and clarity, we also define. 



^eff = ^ X! X! X! S[n,Mf,Ma] Wf, Pa, L) (7) 

N Mp Ma 

Here N is the number of the closed loops in the temporal direction {e.g, 
in Fig. 1 N=2, 4 and 6 respectively). For SU{2), it is an even integer imply- 
ing the invariance under the Z2 center symmetry. Si^^MfMa] Wf^ (^a, L) is a 
functional of the Polyakov loop L{n). Mp and M4 are integers and denote 
the highest powers of (3p and (3a in S\^m,Mp,MaM- interested in 

investigating the presence of a tricritical point in the above theory, we com- 
pute the polynomial S\^N^MF,MA]iPF-, Pa) for N=2,4,6 in the strong coupling 
expansion. The higher order terms in L are found to be negligibly small in 
the region of interest in the extended coupling plane and do not play any 
significant role in determining the phase boundary and the nature of the 
deconfining transition and therefore will be ignored here onwards. 

t A generic gauge invariant term constructed out of Uo{nQ,n) at spatial site (ft) is 
Tr (W^ (tt-))"^ for some integer m. After making a gauge transformation to diagonalise 
W^iji), the above term can be written in terms of L{n). 
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As mentioned earlier, the measure over the Polyakov loop can be fixed by 
the gauge invariance and periodic boundary conditions at finite temerature. 
We define, 



U^in) = U^in)cro + iY,Ul{n)a, (8) 



i=l 



In the above equation, (Tq and CTj are identity and the Pauli matrices 
respectively. Now at a particular spatial site (n), the measure over temporal 
links, ni^'=i / dUo{no, n), can be written as n^'=2 / dUo{no, n) J dU^{n). As 
the integrand at site (n) depends only on UQ^n), the first (A^,- — 1) Haar 
integrals can be trivially done and the integrals over U[{n) give the Jacobian 
Jin) 



Jin) = Q{l-L'in)y (9) 

Here the factor ^ in the measure is due to the normalisation (/ dU = 1). 

The leading strong coupling diagrams contributing to the effective action 
are shown in Fig. l-a,b,c. We define = ^ ^ ^^^^ some 

extensive computations ^ for Nr > 2, the leading local contributions to the 
effective action are: 



^[0,4,2] = (1F + IA + iFlA - ^ [ipf + ^ [lA? j Np (10) 



4' 



=4(^1 Y.L{n)L{n + ^) (11) 



n,i 



S[4,2N.,N.] = ilF + IaT^ Y {L{n),L{n + z)) (12) 

+1^aPf\ Z{L{n),L{n + t)) (13) 



^The details will be presented elsewhere. 
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In the above equations we have used the notation: 



Y{L{n),L{n + i)) 
Z{L{n),L{n + i)) 



(^) Y.{^^L\n)L\n + t)-%L{n)+A^ (14) 

1 s JVr / 

-j ^ (32L^(n)L3(n + i) + 8L(n)L(n + i) 
(l-2(L2(n) + L2(n + z))) +4] (15) 



The local higher order corrections to the above leading effective action 
come from filling the diagrams in Fig. 1. by extra plaquettes. The two and 
four plaquettes corrections to and terms|§] are, 



'[2,7Vt+4,2] 



^1 



[A,2Nr+2,Nr + l] 



1 (iV. + 2) 

37. + 7A-^- 



+ (iV, + 2)(7^)2 

+ ^iV,(7A)')A^.%iv.,o] (16) 

{if + lAi""-'^ ({if + 7a)' - {if?) 

Y{L{n),L{n + i)) (17) 



^1 



[4,27Vr+4,Ar^+2] 



+ 



(i) ""^^2 ((7F + lA? - (^) (7f)')' (7F + 7a) 



Arr-2 



+ j ^^C'l (^(7F + 7a)' - {ifY - 27^7Aj (7f + 7a 



14\ 



\Nr-l 



Nr[l + \n^) [{if + 7a)' - ) 7I) {IF + 7a) 



^) {Nr + 1) (7F + 7a) 



iVT+2 



F(L(n),L(n + 2)) 



Sin /3yii = and the plaquettes correction term in the effective action is computed 
upto order 7^?, i.e, the first term in (16). However, its coeffecient is "fp instead of — ^7_f 
in our case. 
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For iV > 4 in the strong coupling diagrams also give non-local contri- 
butions to S. However, these contributions are proportional to much higher 
powers of the fundamental and adjoint couplings and are small near the 
investigated phase boundaries and thus have been ignored. The Landau 
Ginzburg effective potential can be trivially obtained by putting L{n) = L 
= constant in the above effective action. The effective potentials for /Sa = 
0.5,0.75,0.9,1.4,1.5 and 1.75 at A^^ = 4 are plotted in Fig. 2-a,b,c,d,e,f. In 
Fig. 2-a,b,c corresponding to a second order deconfining transition, the three 
effective potential curves are at f]'^^^^^"-^ and /3|r***™'±0.05. In the Fig. 2-d,e,f, 
where the transition becomes first order, the potentials are at ^ir***"^"' and 
±0.008. They clearly show a dramatic change in the physical proper- 
ties of the theory for large values of the adjoint couphngs. Both in the second 
and first order regions the sharpness of the transition increases with increas- 
ing values of the adjoint coupling as found earlier by Monte Carlo simulations 
[ pO[| . Moreover, in the first order transition region, the discontinuity in the 
Polyakov loop is large, leading to strong first order transitions. In fact, at 
/3a = 1-5 the Monte Carlo value of < \L\ > is 0.5 pO], in close agreement with 



its analytical value 0.55 from Fig. 2-e. The predicted for the various 

values of the adjoint couplings Pa at A^t = 4 along with their corresponding 
Monte-Carlo values [^, ^ are given in Table 1. In all cases good agreement 
is found between the two. The tricritical points for A^^ = 4 and Nr = 6 
are found to be at [/3f = 1.352, /3a = 1.296] and [pp = 1.6711, /?a = 1-864], 
respectively. For Nr = 2, the above contributions to the effective potential 
are slightly different and are presently under investigation. However, both 
Monte Carlo simulations and lower order effective action indicate that the 
tricritical point is in the vicinity of Pa ~ 0.7 — 0.8. The above strong coupling 
shifts in the direction (Pf = oo, Pa = oo) is probably an indication that the 
continuum limit is not affected by the existence of the tricritical point. How- 
ever, such a contrained movement of the tricritical point with N^-, if also true 
in the intermediate and weak coupling regions, requires analytical explana- 
tion. The more interesting possibility of the existence of a new continuum 
theory in the extended coupling plane is still an open question. It will also 
be interesting to investigate other quantitative properties (^e.g, of this 
model for higher values of Pa by Monte Carlo simulations |2^ . This requires 
very large lattices because of the known strong violation of the scaling rela- 



tion in the large Pa region E^, ITSf. However, the more important problem is 
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to understand the origin for such a drastic change in the quahtative proper- 
ties of the theory on the lattice itself and to identify the degrees of freedom 
and the mechanism which make the deconfining transition first order. An 
obvious guess for the above phenomenon is the different global properties of 
the SU(2) and S0{3) groups. This difference can be formally described by 
writing the 5*0(3) part of the action (|^) in its Villain form [|1^. In the latter 
model the above difference corresponds to the Z2 vortices associated with 
plaquette fields, ap = ±1 in the action. In the extreme case, when all the 
vortices are suppressed the above Villain action reduces to the Wilson action 
with modified coupling with only a second order transition. Therefore, it is 
expected that the above change in the qualitative properties is due to the 
condensation of the vortices. A more careful study both via strong coupling 
as well as Monte-Carlo simulations is possible by controlling these objects 
with a potential term A J2piaquettes c"p the action and studying its effect on 
the deconfinement transition. It will be interesting to see the dependence of 
the tricritical point on A. Work in this direction is in progress and will be 
reported elsewhere. 

Table 2 gives the predicted critical couplings and their Monte Carlo values 
^] along the Wilson line {(3a = 0) at different A^,-- It clearly shows that the 
predictions of the strong coupling analysis for the deconfining transition are 
very close to the corresponding Monte Carlo values up to A,- = 16. Further 
comparisons in the region even closer to the asymptotic scaling region could 
not be made because of the lack of Monte Carlo data. The above validity of 
the strong coupling expansion in the scaling region (A^^ > 4) found by Monte- 
Carlo simulations and beyond (Sp = 2.2, where the crossover from strong 
coupling to weak coupling is expected, is again a surprising phenomenon. 
However, it needs further verifications with other lattice observables and a 
careful study of the crossover region at T 7^ 0. 

In the forementioned recent paper ||25| studying SU{3) deconfinement 
transition at non-zero temperature, no evidence for the change in the order of 
the deconfining transition was reported. Also by extrapolating the Polyakov 
loop data to infinite volume, the deconfinement transition was seperated 
from the discontinuity in the plaquette corresponding to the bulk transition. 
However, scaling laws for the plaquette discontinuity remained untested. The 
discrepancies in the qualitative aspects of the phase diagrams of SU(2) and 
SU(3) lattice gauge theories need an explanation. One should again check the 
role of Z{3) vortices in SU{3) deconfinement transition. In the SU{2) 
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plaquette susceptibility peak was found to decrease on going to relatively 
large volume. Therefore, it is important to to verify the scaling behaviour of 
the bulk transition in the case of SU{3) also. 

Finally we summarize the findings of this paper. We have analytically 
confirmed the earlier Moonte Carlo evidence for the existence of the tricritical 
point in the extended SU(2) model. All the qualitative as well as quantita- 
tive features of the deconfinement transition in the extended coupling plane 
found by Monte Carlo simulations are reproduced to a good accuracy. As 
discussed above, this study also opens some new analytical as well as nu- 
merical avenues for understanding and analysing the new features associated 
with the extended model. Surprisingly, it is found that the strong coupling 
predictions for the critical couplings are valid well within the scaling region. 
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Table 1 



Table 1: The values of the critical couplings (A^,- = 4) at various (3a- 



Pa 


j^criticai (^Qfj-Qj^g CoupUng) 


^critical i^Monte Car/ot^l) 


0.5 


1.793 


1.83 


0.75 


1.619 


1.610 


0.9 


1.531 


1.489 


1.1 


1.432 


1.327 


1.4 


1.314 




1.5 


1.28 


1.05 


1.75 


1.147 





Table 2 



Table 2: The scaling of the critical coupling with A^,- at (Sa = 0. 



Nt 


Identical (^gi^Qjig Coupling) 


f^criticai ^Monte Carlo^'^^) 


4 


2.00 


2.2986(6) 


5 


2.25766 


2.3726(45) 


6 


2.38535 


2.4265(30) 


8 


2.51942 


2.5115(40) 


16 


2.81613 


2.7395(100) 



[1] The data taken from ^ . 
[2] The data taken from 
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1 



= 2 



Figure 1: The leading SU(2) strong coupling diagrams contributing to Z^, L"^ and 
terms in the effective action, respectively. 



16 





-0.6 -0.4 -0.2 0.2 0.4 




-0.6 -0.4 -0.2 0.2 0.4 




Figure 2: The effective potentials at different adjoint couplings and = 4. The 
figures [a] = 0.5, [b] /3^= 0.75, [c] = 0.9 and [d] = 1.4, [e] = 
1.5, [f] f^A =1.75, correspond to second order and first order deconfinement 
transitions, respectively. The three second order (first order) effective 
potentials in each figure are at ± 0.05 (0.008). 
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